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C_) ' Abstract 

We present a rigorous and rather self-contained analysis of the Verdet constant in graphene- 
like materials. We apply the magnetic perturbation theory to a tight-binding model and obtain 
£SJ ■ a relatively simple and exactly computable formula for the Verdet constant, at all tempera- 

tures and frequencies. Moreover, for our nearest-neighbor model of graphene we show that 
the transverse component of the conductivity tensor has an asymptotic Taylor expansion in 
the external magnetic field where all the coefficients of even powers are zero. 

Si 

1 Introduction 

Faraday rotation is a dispersion effect, discovered around 1845 by Michael Faraday. It consists of 
the rotation of the polarization plane of a linearly polarized light-beam passing through a material 
■ in the direction of the applied magnetic field. The Faraday rotation angle $ is defined as 



>< 



_ u)d(r]- - 7? 4 



VO ■ 2c 
(N , 

where d is the thickness of the material, while r\- and 77+ are respectively the refraction indices of 
the right and left circularly polarized radiation of frequency oj. As we will explain in (jl.2l) . if the 
applied magnetic field b is small, then the rotation angle can be expressed as $ = dbV where V is 
the Verdet constant, named after Emile Verdet pQ who was among the first ones (around 1865) to 
advocate the use of Maxwell's equations in explaining dispersion. At its turn, the Verdet constant 
depends on how the transverse conductivity coefficient behaves as a function of b near 6 = 0. 

In this paper we put to work the general method proposed in [2] and apply it to the case of a 
tight-binding model. One advantage of a discrete model is that we can give much shorter and less 
technical proofs for both the thermodynamic and adiabatic limits. Another advantage is that the 
structure of the Bloch bands is much simpler in the discrete case, where in many interesting cases 
we know the exact expression of the fiber Hamiltonians, given by finite dimensional matrices. 

1.1 A description of the Faraday effect 

Let us briefly discuss the physical problem following [3J and g] . Starting from the classical Maxwell 
equations, one can derive the Faraday rotation angle for a two-dimensional quasi-free electron gas 
placed in a magnetic field B (of strength 6) applied perpendicular to the layer. The incident 
light-beam is monocromatic. The electrons are confined within a slab of thickness d and move 

freely along the xy plane. Denoting by | cr jfc D ''| . t ne three-dimensional conductivity tensor of 

the slab, and by Oj_ 3D ' = er^ D ' ± ia^\ the complex refraction indexes (77 — in) for the right (+) 
and left (— ) circularly polarized light are given by 

(V± ~ iK ±f = M e f 1 ~ (47ri/we)c4 3D) l . (1.1) 
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The element cr^ 3 ^ does not contribute to the Faraday rotation. Considering a finite thickness d of 
the slab, the three-dimensional conductivity tensor is related to the two-dimensional conductivity 
tensor by the following expression 

-r=-r^ j,* €{i,2}. 

Let us define the real- valued function 



Air \ 2 /4tt n 2 

1 H u + — v 

we / \ ue 



/ 1 

cos — arctan 
I 2 



where |u| < j^,v € R. Now we can express rj± appearing in expression (|1.1[) (to be chosen 
non-negative from physical considerations) as: 

/ — f f-U ( 2D ) ( 2D ) 

»7± = VMe/ (±^2i 



We are interested in studying how the Faraday angle $ behaves as a function of the strength of 
the external magnetic field, b. Let us assume (we will prove it later) that one can write down an 
asymptotic expansion of <j^" (b) in powers of b. Moreover we assume that for the off-diagonal 
tensor element we have: 



while 4i D) (&) can be expressed as: 
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Then 



0(6) = -d6 



c du 



(O,ol¥)+O(b*). 



(1.2) 



Thus the Faraday rotation angle is linear in b near zero, and the linear term can be put in the 
form i? = dbV, where V is the Verdet constant. Calculating the Verdet constant is therefore a 
question of finding a computable formula for the coefficient in the asymptotic expansion of 
the off-diagonal conductivity element cr^K 

The structure of our paper is as follows. In the rest of this section we introduce the mathemat- 
ical notions which are necessary in order to properly define the transverse conductivity element in 
(I1.15[) . In Section [5] we formulate the main result of the paper, namely Theorem 12.11 Sections [31 
S] and [S] contain the proofs of the three statements of our main theorem. In Section [5] we present 
our conclusions. 



1.2 The configuration space 

We neglect the electron-electron interactions. To simplify notation we work in a system of units 
where h = 2rn c i cc t r0 n = e = 1. Let a\ = (ai,0) and 0,2 = (0,02) be two orthogonal vectors in R 2 . 
Define the Bravais lattice, 

r = {7 e R 2 : 7 = moi + na,2, m, n G Z}. 

The basis f2 is modeled by i/q points whose position vectors are denoted by {^ n }^=i; one of them 
is always the origin, thus also belongs to T. For graphene = 4 if we want to have a 'straight' 
Bravais lattice (see figure H} . 

The total configuration space is denoted by A and describes the ion-core mesh. A point in A 
is called a site. Moreover, we have the identity: 

A = r + n = {~f + xeR 2 ■. 7 e r, x e fi}. 
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Define the reciprocal primitive vectors {b n }^=i as: 

a, • bj = 2ir5ij, 

where Sij is the Kronecker delta. The reciprocal lattice (dual lattice) T* is defined by 

T* = {7* e R 2 : 7* = m&! + nb 2 , meZ.neZ} 
We denote by fT the first Brillouin zone for the dual lattice: 



n* = < tibi +t 2 b 2 



1 1 

- < U < - 

2 ~ 2 



The one-electron Hilbert space is £ 2 (A). In the absence of the magnetic field, each electron will 
be described by a one-particle Hamiltonian Ho, where Ho is a nearest-neighbour tight-binding 
operator. If we denote by {5x}xeA the canonical basis of £ 2 (A), then H has a kernel 



h (x,y) = (S x ,HoSy) 



(1.3) 



which is zero if \\x — y\ \ is larger than some constant. 

Graphene is one material which can be described in this way. Since it is important to have a 
'straight' Bravais lattice, we need to extend the minimal model to a basis with four sites instead 
of two. Let a > 0. We let the basis consist of four atoms Hq ■= placed at positions 



(0,0), 



,(2) 



(a,0), x 



(3) 



3a y/3a 



„(4) 



5a y/3a 
Y' ~T~ 



(1.4) 



as indicated in figure [TJ The Bravais lattice associated with this model is generated by the two 
© - (A) ~ 0) (0{V3a) 





Fi gure 1: A nearest-neighbour model of Graphene. To each unit cell belongs a basis of four sites numbered 
one through four, the grey circles indicate the basis positioned at origo. The full lines indicates the non-zero 
hopping matrix elements ho(x,y), where at least on e of x or y belongs to the origo unit cell. We denote the 
nearest-neighbour distance by a. 



vectors a\ — (3a, 0) and a 2 — (0, V3a), see figured] , and the first Brillouin zone is thus 



n* G = 



IT IT 

3a ' 3a 



\/Sa ' \/3a. 

In this model a site have nearest neighbours either at the three relative positions 

' a \/ia \ I a \/ia \ 



(a,0), 



2' 2 



and 



2' 



(1.5) 



(1.6) 



1 INTRODUCTION 



4 



or at the three relative positions (-a, 0), (§, ^f 2 ^ and , -^f 2 ) . 

This leads to a kernel hg(x, y) which equals (here 5 means Kronecker delta): 



£ S( X1 - yi , (-!)*§) £>(a: a - V 2, (-1)^ J + " (-1)^)^2,2/2). (1.7) 



-\/3a 

fc=i j=i \ / ;=i 

In other words, (x,y) = 1 if ||cc — y|| = a, and h$(x,y) = otherwise. 
If we restrict ourselves to a finite crystal, it will be modeled by 

A N := {(j + x) E A : 7 = mai+na 2 , |m| < N, \n\ < N, x G CI}, N > 1. 

This central region contains | An | = (2N + l) 2 unit cells. The characteristic function of the central 
region is denoted by xn- The Hamilton operator subject to Dirichlet boundary conditions (DBC) 
in Ajy is: 

Ho,n = XnHoXn- (1-8) 

In general, if some operator O initially defined on the whole site-mesh A is afterwards restricted 
to A^v, we denote this restrition by On = xnOxn- 

Now we include a constant, static and external magnetic field into the model, which is thought 
of as having always existed (before the light pertubation is turned on). The magnetic field is 
assumed to be perpendicular to the layer, directed in the z+ direction, having constant magnitude 
b. All our vectors can be seen as three dimensional, and we associate with the symbol x both the 
2d vector (xx, X2) and the 3d vector (xi, X2, 0). This enables us to use the cross-product shorthand 
x x y for 2d vectors. 

Invoking Peierls substitution, our magnetic Hamiltonian Hb in a constant magnetic field has 
an integral kernel hh{x,y) which is obtained by multiplying ho(x,y) by a phase factor 0[6j[7], 

h h (x,y) = (6 X7 U h S y ) = e ib ^ x ^ho(x,y), 

with ^ 

v{x>y) = -^(yiX2 - xm) = ^ [G/i>2fc>°) x ( a; i> a; 2,o)] z , (1.9) 

where [v] denotes the z-component of the 3D vector v. Note that ip is anti-symmetric, that is 
ip(x,y) — —ip(y,x). We denote by fl(cc, x' , x") the magnetic flux of a field of unit field-strength 
through the triangle generated by the sites x, x' and x" . 



&(x,x',x") = 



\[{x'-x")x{x-x') 



We note that the following identity holds: 

fl(cc, x', x") = ip(x, x') + ip{x', x") + cp(x", x). (1.10) 

Finally, we define the position operators Xi : £ 2 (A) £ 2 (A) and X 2 : £ 2 (A) —> £ 2 (A) by their 
action on the basis elements: 

X„6 X := x v 5 x , Vx = (xi,x 2 ) e A, ^e{l,2}. 

1.3 Derivation of the conductivity tensor from the Kubo formalism 

We use the same strategy as in [8] and [2] in order to express the conductivity tensor as a local 
trace. 

We denote by /j, the chemical potential and by (3 = -^Kp the inverse temperature. In the remote 
past our system is at thermal equilibrium, described by the Fermi-Dirac density operator: 

f?0,N = ./FD(H b ,N), fFu(z) 



£(z-M) + 1 ' 
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The system is perturbed by an incident monochromatic light-beam with a complex frequency 

lo: 

to = loq — iji, Re (a;) = loq > 0, Im(w) = — rj < 0, (1-11) 
and whose electric field is parallel with the x-axis and given by 

E(t) = {E x (t),E y {t)) = (Re (Ee lult ) , 0) = Ee 71t (cos(w <), 0) , E > 0, t < 0. 

This electric field generates an external time-dependent potential term V B (t)- The full time- 
dependent Hamilton operator is now given by 

H B (t) =H b +V E (t), V E (t)=Re(Ee^ t )X 1 . (1.12) 

To find the density operator at an arbitrary time t < 0, we need to solve the dynamic (Liouville) 
equation with the initial condition at to = ~~ o° : 

Ll ' N (t) = [H EiN ,e BN (<;)], lim e EN (t) = / F D(H b ,N)- (1-13) 



-oo 



Note that at finite Af all operators are finite dimensional matrices. The negative imaginary part 
of ui has the effect that Ve(-) is absolutely integrable near t = — oo. Therefore equation (|1.13l) 
has a unique solution which can be iterated. Now we will identify the linear response coefficients. 

1.3.1 Current density 

The current operator in direction v is defined as 

ji/,b,N '■— *[He : n , X„ jN ] = i[H b ,N , X Wj n] v = 1, 2, 

which is independent of -E because V B (-) commutes with X„.n- The current density in the y- 
direction at t — is given by 

J2,b,N(E) = ^Tr{p E . N (t = 0)j 2 , b , N }. (1.14) 

\ A N\ 

It is well known that J2,b,N(E) has an analytic expansion in E near zero, and we will show that 

J2,bA E ) = E ^i(b, N) + 0{E 2 ). (1.15) 

The linear term in E in formula (|1 . 15|) defines the off-diagonal conduction element that we seek 
to identify. 

2 Main results 

The operator norm of Hb,N is bounded from above by a constant uniform in b and N, thus its 
spectrum lies in a sufficiently large closed interval I, uniformly in b and N. 

We define a smooth, simple and positively oriented closed path which encloses the above 
interval while oo lies outside ff. Moreover it has to be close enough to the real line such that /fd(z) 
has no singularities inside % see fig. [2] 

As a preliminary result, we shall see that a 21 can be put in the form: 

V_ 

(v 2 +^ )\An\ 



*2i(b,N) = ~ ,„ a , . a VA . , Tr{i[j 2|b|N ,X 1 , N ]/ FD (H N , b )} (2.1) 



+ Re 7; — 11. — 1 / dz /fd(^) (Tr{(H Nib - z + u) 1 ji,b,N(H N ,b - z) l h,b,n\ 
2ttuo\K n \ Jcg V I J 



+ [z Z + CO , 



where [z — > z + 00] denotes a trace like the preceding one, just with z substituted by z + to. For 
z € p(H b ) define the operators on £ 2 (A): 

D b (z) := (H b - 2 + w) _1 ji.b(H b - «) _1 j2,b (2.2) 
where j^b := i[Hb, X„], v = 1,2. We are now ready to state the main result of this paper. 
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Figure 2: The path V should enclose I, and u> should be outside ' 



Theorem 2.1. Assume that the real part of the frequency loq is large enough. Then the following 
statements hold true: 

(i). Assume that c €§ contains the spectrum ofHb, while ± ujq does not. Then the transverse 
conductivity admits both the thermodynamic and the adiabatic limit, and we have: 



o~2i(b) := lim lim er 21 (&, N) 

i)->0 JV->oo 



Re ^T / dz /fd(^) Yl {B h (x,x;z) + T> h (x,x;z + uj )) 



(2.3) 



(ii). In our model of graphene (see (|1.7|) ). 02^(0) = 0. 

(Hi). The function b t— > <J2i(6) is smooth and has an asymptotic expansion in b around 0. All 
the derivatives of o~2i at zero can be written only in terms of the fiber operators associated to the 
Bloch decomposition of Ho . In particular, for our model of graphene, all even Taylor coefficients 
are zero: 

d n ag 



-(0)=0, n = 2p, peU. 



(2.4) 



3 Proof of Theorem l2Jli) 

3.1 Derivation of formula (12. ip 

If A is a bounded operator, then we denote its expression in the interaction picture with 

A(t) = cxp (iiHb.N) A exp (— jiHb,N) ■ 
By standard perturbation theory, we can write 



£»n,e(0) = /FD(H N , b ) - i 



V n ,e(s), /F D (H Nib ) +0(E 2 ). 



Thus in the linear response approximation we set 

1° r- 
Qn e nn(0) = /FD(H N , b ) - iE / dsRe (e^ s ) X 1N , E (s) , / FD (H N , b ) 



(3.1) 



The value of the current density in the y-direction in the linear response regime is given by 
J 2 ,b,Nji„(£, * = 0) = Tr { ( / PD (H Nib ) -iE f ds Re (e ius ) [x 1N)E (s) , / FD (H N .b) 



j2,b,N 



N 
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Using the trace-cyclicity rule Tr{[A, B]C} = Tr{B[C, A]}, the equilibrium current J2(0) is shown 
to equal zero: 

J 2 ,kn(0) - Tr{[/ FD (H Nib ),H N , b ]X 2)N } = 0. 

By examining the formulae (|1.14p . (|1.15p and (|3 . 1 1) we can single out the transverse conductivity 
term: 

1 r° r~ i 

o-2i lN (i = 0) = -j— / dsRe (e lu]s ) Tr{i X lib , N (s) , / FD (H N . b ) j 2 , b , N }. (3.2) 

\^N\ J-oo 1 J 

The trace in the integrand of equation (|3.2[) is a real number, thus formula (|3.2|) can be re-expressed 
as 



C21,n(£ = 0) 



-^Re / dse^ s Trjz^.^N, e' ;sHN . b X liN e- 4sH ^]/ FD (H N!0 )}. (3.3) 

l A AM J-oo 



By partial integration we re-express (|3.3p as 



We note that 



( V 2 + ^ )\A N \ 



Tr {i[j 2 ,b,N, Xi iN ]/ FD (H N!b )} 



(3.4) 



Re Tiri I d S -e^-^(Tr{*[j 2 , b ,N,e lsHN ' b X liN e- 4sH ^ b ]/ FD (H N . b )}). 



|A 

d_ 

df 



Ji,b,Ne 



by the definition of j ljb . By using the trace-cyclicity and by noticing that / F D(HN, b ) and 
e is(w+H N b ) commute, we have 



o-2i,n(£ = 0) 



■ Tr {i[j2,b,N, X liN ]/ FD (H Nib )} 



+ 



(V 2 +^ )\A N \ 

H^Rc / f-Trf e - HNb )j l!b , N e- 4S ( HNb -"V FD (H Nib )j2, b ,N} 

\ A n\ J-oo \0J ^ > 

- Tr {e 4s(w+HN ' b) /FD(H N , b )j l!b , N e 4sHN b j 2 ^ N )}) ds. 

We can use the Cauchy integral formula to express the operator e ,s ("+ H N,b)j' FD (ij N|b ) by a curve 
integral in the complex plane, involving the resolvent of HjM, b : 

e «("+HN, b ) /FD (H N!b ) = ± I e -(-+")/ FD (z)(H N , b -z)- 1 rf 2 , 

where the path ^ encloses, but has no points in common with, the (real, bounded) spectrum, 
(j(H b .N)- We claim that it is possible, given and fi, to choose such a curve such that w 
lies outside ff, and so "close" (small imaginary parts of z e ^) to the real line, that /fu(z) = 
(oPiz-p-) _|_ h as no singularities inside ^ . See figure [2] This leads to 



0"21,n(* = 0) 



(V 2 +^ )\A N \ 
( 



Tr{i|j2 lbl N,Xi iN ]/ FD (H N!b )} 



1 



+ Re / ds- 

-oo U\A N \ 



sHN ' b ji,b,N (± £ dz e~ is ( 2 -^/ FD (z)(H N , b - z)- 1 ) j 2 , b , T 



"(H Nib -^) /FD(HNb) 



Tr < 



I 

2^ 



d2e fa (* +w ^ D («)(H Nib -«)- 1 )j 1 



,b,N e 



' j2,b,N 



-„"( H N,b + ") 



/FD(H^ 



(3.5) 
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The two integrals, §^ dz- ■ ■ and J_ ao ds- ■ ■ in equation Q3.5P are both absolutely convergent, there- 
fore we can exchange integration order (Fubini). Furthermore, as e ±ts( - z±u] is nothing but a 
complex scalar, we can freely place this factor in the operator product, 

CT21,N(f = 0) = ~ 7 9 , ?m T Tr {i|j2,b,N, Xl n]/fd(Hn, b)} 

(ri 2 + w£)\A N \ 

+ Re / dz / d,s— ^ ( Tr (e"( HN -*-* +u >j libJ ,;^/ FD (z)(H N . b - z)~% h A 
- Tr |^-/ FD (z)(H N , b - z)- 1 Ji,b,Ne Js( ^- H - b) j 2 ,b,N|) . 
Integrating with respect to s we obtain (formula (|2.1|Q : 

<T 21,n(* = 0) = — , „ 2TH 1 ^ {*[j2,b.N, Xi N ]/ F D(H N ,b)} 

+ Re/ dz o j , /FD(^)(Tr{(H N .b-z + ^)" 1 j 1 .b.N(H N , b -2)- 1 j 2 ,b,N| (3-6) 
+ Tr {(H N;b - z) _1 ji,b,N (H N , b - z - j2,b,NM • 

3.2 Off-diagonal localization for resolvents 

In this subsection we only work with operators defined on £ 2 (A). It should be understood that 
the results (|3.1I) - (I3.5[) also hold true, even if A is replaced with An, uniformly in N. 

Definition 3.1 (Schur-Holmgren bound). For a linear operator A E B (£ 2 (A)) with a kernel 
a(x,x'), we define the Schur-Holmgren bound by: 



:= max < sup }^ \a(x,x')\, sup y^\g(x,x?)\ \ 



(3.7) 



If an operator has HA^ < oo it is said to be Schur-Holmgren bounded. 

The following bound is well known, and we give it without proof: 

Lemma 3.2. Let A € B(t 2 (A)). If \\A\\ 1 < oo, then ||A|| < ||A|| 17 where ||A|| is the usual 
operator norm. It holds that 

l|A|| < ||A|| r (3.8) 

Definition 3.3 (Exponentially Almost Diagonal Operator). Let A e B (^ 2 (A)). We say that A 
is exponentially almost diagonal, if there exist two constants C\,Ci, both strictly positive, so that 
the kernel of A satisfies 

\a(x,y)\<C ie - c ^ x -ytt (3.9) 

for all x, y in A. 

The proof of the next lemma is straightforward and thus omitted. 

Lemma 3.4. An exponentially almost diagonal operator is Schur-Holmgren bounded (and by 
lemma [37B bounded). 

We now show a property for exponentially almost diagonal operators, which is a much simpler 
version of the Combes- Thomas estimate for resolvents of continuous Schrodinger operators [§1 HH] ■ 

Proposition 3.5 (CT-property). Let z £ p{A) where A £ B (£ 2 (A)) is a self-adjoint exponentially 
almost diagonal operator. Let constants G\ and Gi be defined as in definition (|3.9p . 
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Then the resolvent (A — z) 1 is also exponentially almost diagonal. That is, there exists two 
positive constants C3 and C4 such that the kernel of (A — z)^ 1 fulfils 

\(A-z)-\x,y)\<C 3 e- c ^ x -yK 

Consider the situation where z is restricted to a closed subset ^ C p{A). Then C3 and C4 can be 
chosen uniformly in dist( c #', cr(A) ). 

Proof. We only sketch the main ideas, see [IT] for a related result. For a > 0, and a fixed lattice 
point x a £ A, define the operator A a : £ 2 (A) — > £ 2 (A) by 

A Q := e a ll-- a; o|l Ae -a||--a;o||_ 

It can be shown that 

a -> 

HA-AJ! -0. (3.10) 

By the identity 

(A a - z) = (I - (A - A a )(A - «)- x )(A - z), 
the operator (A Q — z) is invertible for a small enough depending on z. We can write the equality 

e -a||— sell ( Aa — z) = (A — z)e- Q "' - X °H , 

which implies that 

(A - z)- l e-"^ - x "^ = e— "H- x (Aq _ 

and which gives 

e a\\- -X \\ ( A _ z J-l e -a||- -Xo|| = ( Aq _ 

If we apply both sides on the basis element 5x and then take the scalar product with some Sx 
we get: 

(5 X , e a|1 ' (A - ^-^-"H-^ll^) = e «ll*-*o|| (A _ z )-i( Xj Xq) = {Sx , ( Aq _ z)- 1 ^), 
which means that 

go||a;-»oll I ( A _ «)-i( XjXo )| < ||(A Q - z)- 1 !) . 

This formula can be used in order to get the desired estimate for the integral kernel. We do not 
give further details. □ 



The last preparatory results are concerned with the magnetic translations. We define the 
magnetic translation operator as the operator that transforms xj) £ £ 2 (A) according to the rule: 

(T b , 7 -0) (a:) := e ib ^ x ^iP{x - 7), for all 7 e T, x G A. 

The magnetic translations have certain properties which we list here without proof. The inverse 
of the magnetic translation operator obeys 

7b, 7 — 7b, -7- 

The Hamilton operator Hb commutes with Tb.-y for all 7 € T. For z G p(Hb) we have 

(H b - z)- 1 = Tb, 7 (H b - z)- Y T h .^. 
The same property is true for the current operator whose integral kernel is 

j„,b(x, y) = i&bfo y)(i/„ - x v ) = e ib ^ x Mj v>0 (x, y), v e {1, 2}. 
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3.3 Proof of the thermodynamic and adiabatic limits 

Let us briefly discuss the first term appearing in (|2.ip . that is 

12, ~2 \ I A 1 Tl '{*[j2.b,N,Xi N ]/fd(H N b )} • 

We have the identity: 

«[i2,b,N,Xi jN ](x,y) = ixN(x)h b (x,y)(y 2 - x 2 )(yi - xi)xn(i/), 

which defines a bounded operator due to the localization properties of ho. This means in particular 
that Trji^bjNiX^Nj/FDlHi^tOIVI-A-ivl * s bounded in N uniformly in r], thus after the adiabatic 
limit this term will disappear anyway. That is why we only treat in detail the second term of 

First, we need to introduce some notation. Let < e < 1. We divide our finite box, Ajv, 
into an edge region A N of width [N £ ] unit cells, and a remaining core part, A N , see figure [31 We 
have that An = A^v \ A^. For practical reasons we work with N e instead of its integer part. The 




2N+1 



Figure 3: Ajv is split into a core part, Ajv, 2(N — N € ) + 1 unit cells wide and an edge region Ajv of width [N*] 
unit cells. < e < 1 



number of sites in An is proportional with N 1+e , which means -j$r ~ * as N tends to infinity. 
As a consequence, the number of sites in An behaves like AN 2 for large N. 

3.3.1 Geometric perturbation theory 

To simplify notation, we introduce a shorthand for the characteristic functions xn = Xajv > 
Yd = y= . We now introduce an auxilliary operator by: 

Ab,N(z) = Xn (Hb - z)' 1 xn + Xn (Hb,N - ^y 1 xn, z € c i . 
If we multiply A^n^) on the left by (Hb,N ~ z ), we have 

(H b ,N - z) A b ,N(^z) = (H b ,N - z) xn (H b - zy 1 xn + (H bjN - z) xn (H b , N - xn- (3.11) 

Note that for large enough N, the distance between An and A \ An becomes larger than the 
interaction range of H bj N- This implies that (H b .N — z ) xn = (H b — z) xn- Therefore, if we use 
this in the first term of formula (|3.11[) . and using that xn + Xn = Xn, (|3.11j) becomes: 

(Hb,N - z) Ab,N(^) = XN + [H b , xn] (Hb - z)^ 1 xn + [H b:N , Xn] (H bjN - z)~ x xn, 
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which is equivalent with: 



where: 



(H b , N -z) 1 = A blN (z)- (H b , N -z) 'Bt.Nfz), (3.12) 

B b ,N(z) = [H b , xn] (Hb - Xn + [H b , N , xn] (H b ,N - z)' 1 xn- 

We insert (j3. 12[) in formula (j3.6|) obtaining several terms. We claim that only the following term 
contributes in the large N limit: 

Re 2mv\A^~\ ff<# dZ f pD ^ ( Tr {* N *- Hb ~~ Z + UJ V 1 XNji,NXN(H b - z)~ 1 xnJ2,n} + [z -> z + u)]j . 

~"~ ' ' ' ' (3.13) 
The other terms in the expansion of formula (|3.6p have factors of type xn (H b ,N — z) , [H b , xn] 
or [H bj N,%N]- In the large TV limit terms having these factors vanish, which we explain in the 
following. 

To begin with, let us choose one such term which up to trace cyclicity can be written as 

' Tr{[H b ,x N ]Q 6 ,N(z)}, ^e^, 



|Ajv 



where Q;,,n(*z) is bounded uniformly in b, N, and z. Because the operator H b is short-range, we 
can find a projection Pn whose corresponding subspace has a dimension D ~ N 1+c , such that 
PN[H b ,XN] = [Hi»Xn]- We then use the inequality 

|Tr{[H b ,XN]Q b , N (z)}| < CTr{P N }~iV 1+e , 

where C is a constant uniform in 6, N , and z. Thus by dividing by -/V 2 , it will converge to zero. 
The same type of proof applies for all other terms containing xn which can already play the role 
of P N . 

We now consider the term given by formula (|3.13[) . and we want to show that we can replace 
XnJi.nXn with ji. If we write 

Xn =1- (I-Xn), 
then the trace in (13 . 1 3[) can be written as: 

Tr {xn (H b - z + u)- 1 (I - (I - xn)) ji (I - (I - Xn)) (H b - z^XnMn} + [z -> z+u]. (3.14) 

This trace, and thereby formula (|3.13[) . can now be expanded into several terms, all but one 
containing at least one factor of the type (I — Xn)- The terms containing (I — xn) = (I — Xn + Xn) 
as a factor can be written (up to trace cyclicity) in the form 

Re } / dz/ FD (z)Tr{(I-XN)E b (z)xNQi 1 N W+XNQi 2 N (^)}, (3.15) 
2ituj\A n \ Jcg v J 

where E b (z) is an exponentially almost diagonal operator and q[^(z) and q[ 2 ^(z) are bounded, 

~ (2) 

uniformly in b, N and z. The term xnQJ, n( z ) vanishes as N tends to infinity by previous argu- 
ments. Now consider the term containing (I — XN)E b (z)xN- The kernel k(x, y) of this operator is 
zero unless \\x — y\\ > N e , see figure [3] One can easily show that the trace-norm of this operator is 
bounded from above by Cie~ C2N , for some positive constants C\ and C2. Thus this term will not 
give a contribution to the thermodynamic limit. The only remaining contribution from formula 
(13131) is 

Re o — H — 1 / dz /fd(z) (TrjxN (H b - z + w)" 1 ji, b (H b - z)~ 1 j 2 ,bXN} 

27Tw|AAr| J<g \ I J (3.I6) 



Tr|xN(H b -z) x ix,h (Hb - z - u) 1 fabXtiY) 
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Using the operators denned in (|2.2p . the previous formula can be re- written as 

Re^— &&/fdO?) I ^2 T> h (x 7 x;z) + T> h (x , x; z + uj) J . (3.17) 

* \xeA N ) 

Dij(z) is a product of operators which commute with magnetic translations. This implies that the 
diagonal elements of its integral kernel define a periodic function, that is for any x € and 7 G T, 
we have that 

Db,±(x + ~r,x + 7; z) = D bi± (x,x;,z). 
Now the proof of (|2.3p is straigthforward and the thermodynamic and adiabatic limits are proved. 



4 Proof of Theorem I2.1( ii) 

Here we need to compute the transverse conductivity component at b = and prove that it gives 
zero if we work with the operator (|1.7p . We will do this computation for a general nearest neighbor 
model, and use the graphene model only at the end. 



4.1 The Bloch-Floquet representation 

In order to fix notation, we define Jtfp := J®„ d 2 k £ 2 (Vt). The Floquet unitary [TO] U : £ 2 (A) ->■ 
takes vectors ip from £ 2 {A) into J#p, and is given by the well-known formula 

(Uil>)(k,x) = —==^2exp(-ik- ~f)ip{x + j), ken*,xe{l. (4.1) 
V 1^*1 7G r 

If A is any bounded self-adjoint operator commuting with the translations induced by the Bravais 
lattice r, we have that UAU* : 3%? — > is given by 

/■ffi 

UAU* = / a F (fc)dfc, 

Jo* 

where the fibers have the kernels (we assume that a(x + 7, y) has a sufficiently fast decay with 
7): 

a F fey;fe) = ^e- ifc ^a(^ + 7,y). (4.2) 

We denoted the number of sites of fl with vq = \Vt\. For each k £ il* , ap(k) is self-adjoint and 
has a diagonal matrix with un real eigenvalues. We order these eigenvalues in increasing order: 

ei(fc) < e 2 (k) <■■■< e„ n (k). 

In order to have periodic boundary conditions in the fibers, we modify (|4.ip with a complex 
phase and define: 

U F :£ 2 (A)^^ F , (U F ip)(k,x) = J2e~ lk ' { - +1) ^(x + l), (4-3) 

for all tpe£ 2 c (A)). Accordingly U F H U* : Jf F -> J^ F has the fibers 

ho(x,mk) = £ ho(x + 1 .,y)e- lk <^-y). 

If we differentiate the fiber with respect to the first component of k, we have 

^T~ho(x, Vj fe) = - E e~ ik -^-y) l {x 1 + 7l - Ri )h(x + 7, y). (4.4) 

1 7er 
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The expression (14.4)) is nothing but the fiber of the transformed current operator Up (i[H, Xi]) Up. 

In particular, for the graphene-Hamiltonian (jl.7|) we have = 4, with the numbering of the 
basis positions given in (| 1 .4|) : 







2e 2 cos 

2e^ cos 



2 J 









2e^ cos (^*%^ 




2e~ 



We see that all matrix elements are even in &2- The bandstructure is given by the eigenvalues of 
/jo'(fe) and are shown in figure [4] 




Figure 4: The bandstructure of our model of Graphene. Top: a figure of the first Brillouin zone, with the contours 
of e3(fci, &2) plotted in grayscale, also showing some points of interest: G (origo) XI, X2, M, K and K' . K and K' 
are the two "Dirac points" . Bottom: the plot of the Bloch energies along the path indicated in the top figure. 



Now if we have two operators A, B € B(£ 2 (A)), both T-periodic, then it holds that: 

(AB)Gz,y;fc)= £ a(x,x'; k)b(x',y; k). (4.5) 
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When b is equal to zero, formula (|2.3I) becomes: 



°2i(0) = 7777^ Re </ dz /fd(z) ([(Ho-z + cjo) l ji,o(Ho - z) ^2,0] (as, as) + [z -> z + oj }) 

2|0|7ro;o ^ u 

In order to shorter notation, we denote by r (z,fc) — (h (k) — z) -1 the fiber of the resolvent 
(Hq — z) -1 . Using (14. 5[) we can write 



0-21 (0) 



2 1 ^2 1 7TCt^o 



(ro(x,rf\ z - Uo,k)-^ho(x',x!'-, k)r (x", as'"; z, fe)^- fro (as'", as; k)+ (4-6) 

5 £^ 

+r (x , ; z , fc ) — - ft (#' , x" ; fc ) r {x" , a'" ; z + w , fc ) ^- /j (?/" , x; k) 
ok 1 ok 2 

Specializing this formula for h^(k), we see that by differentiating with respect to &2 and then 
integrating with respect to z, the total integrand for the k integral becomes an odd function of 
&2. When we integrate &2 on the symmetric Brillouin zone, the integral giving ct^iW equals zero. 

5 Proof of Theorem 12.11 (iii) 

This section is where the magnetic perturbation theory plays a crucial role. For any z € p(H ), 
we define the operator Sb(z) by its kernel 

s b (x, y; z) = e**™ ((H - z)" 1 ) (x, y). (5.1) 

Notice that a Schur-Holmgren estimate shows that when z is restricted to a compact set ^€ C 
/)(Hb), then ||Sb(z)|| < Co^) for some positive constant Co(^), uniformly in z £ ^ . By denoting 
with I the identity operator, we define 

K 6 (z) = (H b -z)S b (z)-I, (5.2) 

where the operator K^(z) has the integral kernel: 

k h (x, x';z) = e^ x ') Y, (e»* x >V' x '> - l) h (x, y)a (y, x'; z). 

y 

Using the exponential localization of the above integral kernels, together with the estimate 

\&(x,y,x')\ < \\x-y\\ \\y - x'\\, 
then a Schur-Holmgren estimate applied to K b (z) shows that 

sup{||K b (z)||}<6C7(^), (5.3) 

for some positive constant C(f$), uniformly in z G . The constant C( ( ia) only depends on the 
distance between ^ and ct(H ). 

The next lemma is a direct consequence of the above estimates and recovers a well-known 
result about the spectrum stability of H b . We state it here without other details; see [121 [5] for 
much stronger results. 

Lemma 5.1. Let c & C p(H,o) be any compact set. Then there exists b^g > 0, sufficiently small, 
such that c € C p(H b ) for all < b < b^. 
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From now on is the integration contour in the formula giving the conductivity, and z € 
If b is small enough we can write 

(H b - z)- 1 = S b (z)[I + Kb(z)]- 1 = S h (z) - S b (z)[I + Kbfzr^bfz) 

= S h (z)-(H h -z)- 1 K h (z). (5.4) 

We can iterate this and obtain 

(H b - z)- 1 = S b (z) - S b (z)K b (z) + TZ h (z), (5.5) 

with the definition of the remainder term 

TZ h (z) := (Ht-z)- 1 ^. 

Now both factors defining this remainder are exponentially localized, and standard estimates lead 
to: 

sup|7£ b (a;,y;z)| < b 2 c n e~~ c ^ x ~ y ^ , x. y e A, 

for some positive constants en and cn- This shows that the remainder is also exponentially almost 
diagonal, thus there exists a constant Ctz( c ^) > such that 

\\K h (z)\\ <b 2 C n {^). (5.6) 
5.1 The first derivative of 021(6) 

Now we seek to identify the linear part in b of (|2.3[) . Consider again (|2 . 3|) : 

V2i{b) (5.7) 

= Re ^T?n f dz /fd(^) 52 [( Hb ~ z + w o) _1 ji,b(H b - z) _1 j 2 .b + (z ^ z + w )] 

Using formula (|5.5[) . we see that by subtituting (H b — z) -1 with S b (z) — S b (z)K b (z), the error 
we make is of order b 2 and this remainder cannot contribute to the first order derivative at b = 0. 
Therefore: 

bv$ + 0{b 2 ) = Re^^ / dz / FD (z) (5.8) 

52 K S b( z - w o) - S b (z - w )K b (z - w )) ji,b [S b (z) - S b (z)K b (z)] j 2 , b + (z -> z + w )] (x,x). 
xen 

We will now sketch the calculation of the trace over the basis for a given z. The following 
computations hold uniformly in z € *«? and < b < b<g. To simplify notation we write tu instead 
of wo and we introduce the folowing shorthands: 

S := S b (z), S_ := S b (z - w ), "S-type" 
Q := S b (z)K b (z), Q_ := S b (z - w )K b (z - w ), "SK-type". 

For a; € SI consider the element 

[ (5_ - Q_) jib (5 - Q) hb + (z -> 2 + w ) ] (5.9) 

Let us expand it: 

[5-ji 6 5j 2 h - S-jibQhb - Q-hbShb + Q-hbQhb + (z -> z + ui )](x,x). (5.10) 

To show the method of calculation, we first consider the operator product of two of "S-type" 
operators. 

[s b (z - w )ji, b s b (z)j 2 , b ] (as, as), x e n. 
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Here, the 6-dependence of the integral kernel appears only through the exponential phases. De- 
noting by 

(FL)(x, x',x",x"') = ip(x,x') + <p(x',x") + (p(x",x"')+cp(x"',x), 
we see that the above kernel can be written as: 

J2 e «(FL)(x,x' ) x»,x'") Sofe! x >. z _ UQ ) jhQ ( ?B ' ) x ")s {x", x'"; z)j 2 , (x"',x). 

X' ,X" .X 1 " 

It can be easily seen, using (|1.10|) , that 



(FL) = fL(x, x', x") + fl(as, x", x'") 

= 2 [fe-4)(4 -40 - 

The expansion of e lb ^ FL ' > in b is 



\ [{x.2 - 4)(4 - 40 - fei - 4)(4 - 40] + ^ [fe - 40(4 - 4") - - 40(4 - 40] 



gib(FL) = 1 + ^( FL ) + (b 2 (FL) 2 ). (5.11) 

We see that due to the exponential localization of the various kernels, the terms generated by 
0{b 2 (FL) 2 ) will give a contribution of b 2 , thus it can be discarded. The linear contribution from 
the right hand side of formula (|5.11l) is: 

"2~ K—2 ~ X 2)( X 1 ~ X l) ~ fel — x l){ x 2 ~~ x 2)\ + "2"" K—2 — x 2 )( x l ~ X l ) ~ fel — x l)( x 2 — X 2 )] 

^ / / // / / // 

= — [(a; 2 — x 2 )(x 1 — Xj) — (ie-l — x 1 )(x 2 — x 2 )] 

+ "2"" [fe ~ X 2)\ X l ~ x l ) ~ ( x l ~ x l)( x 2 ~ x 2 ) + ( X 2 ~ ^X^l — ^1 ) — (' T l ~ x l)( x 2 ~~ x 2 )] • 

Introduce this into the formula for [<S_ji,b<Sj2,b] (x.x). We now have that the linear term in b of 
[<S_ji,b<Sj2,b] (x,x) is given by: 

ib x - 

2" ^ 

a;',a;",cc"'eA 

[(x.2 ~ x ' 2 ) s o(x, x'; z - uj ) (x[ - x")j 10 (x', x") s (x" , x'"; z)j 20 (x"',x) 

- {x x -x' x )s Q {x, x' ; z - ui ) (x' 2 - x 2 )j 10 (x' , x") s (x" , x'"] z)j 20 (x'",x) 
+ ( x 2 - 4) s ofe x';z- oj q ) j 10 {x',x") (x'{ - x"')s (x", x'": z) j 20 (vt'" , x) 

- (x x - x[)s (x, x';z-ui )ji (x',x") (x 2 - x' 2 )s Q (x" \x'";z) j 20 (vt'" \x) 
+ s (x,x';z- uj ) [x' 2 - x 2 )ji (x\x") (x'{ - x'{') s (x" , x'"; z) j 20 (x'",x) 
-s (x,x';z - u> ) (x[ - x")j w (x',x") (x 2 - x 2 ) s (x" ,x"' ; z) j 20 (x"' ,x)} . 

Switching to fc-space, we have that multiplying an operator-kernel a(x, x') with {x v — x' v ) transfers 
into differentiating the fiber with respect to fc„, v G {1,2}: 

d 

i(x u — x' )a(x, x') — » -7^—CL{x,x[]k). (5-12) 
For example, when computing the local traces we can make the switch 

J2 (■■■)j,Ax,x%--)^ J2 W\ / t d2fc ("0^ofe^;fc)(--0 

•• x,x' -eA ■■■x,x; - en ' ' n * v 

We thus have that the coefficient of the linear term in b of 



Re ^T f dz /fdW V [5_ji,b5j 2 ,b] (x,x) 
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is given by (remember that ro{z, k) is the matrix (ho(k) — z) 



--^r Q (x,x'; z - uj Qi k) -^h (g^ , x" ; k)r (g^' z, k)-^-h (x"' ,x; k) 

+ -^- r o(x,x!] z - ujq, k) 9 ho(x',x"; k)r {x", x'"; z, k)-^-h (x"', x; k) 
ok 1 OK2OK1 ok 2 

- r Q (x, x[ ; z - uj q , k ) h (x' , ; fe ) r (x" , x'" ; z, k) fr ,x;k) 
+ ^-r Q (x,g^; z - cj , k)-^-ho(x' ,x"; k)^-r (x[' ',x'" '; 2, k)-^-h Q (x"', x\ k) 
-ro(x,x';z-L> ,k) 9 h (x' ,x";k) r (x" ,x"';z,k) h (x'" ,x;k) 

OK2OK1 OK 1 Ok 2 

+ro (x, d ; z - ujq , fe) ^ /i (#' , x"; fc) r (x" , x"' ; z, k) h (x'" ,x;k) 
okf ok 2 ok. 



(5.13) 



^2 

Now consider the factor Q = Sb(z)Kb(,z) which is a "SK-type" operator: 
[S b (z)K b (z)] (x,y) 

= J2 e lb ^ x - x 'ls (x, x'; zyW'M (jbfl(x' ,x" ,y) _ A x ») Sq ( x » , y . z y (5.14) 

;c\a;"eA 

We see that e M ( x ,x ^ — 1 is already first order in b, thus we can discard the two terms Q-jibQj2b 
(one from z and one from z — > z + ujq) in (15.101) . when neglecting all terms not linear in b. 
Furthermore, we can reduce (|5.14p to 

e Wx,y) £ So(a . x r. z) [ e m(x>,x»,y) _ ^ ^ x » )sq{x » } y . z)} (5 . 15) 
x',x"eA 

when neglecting higher orders of b. The first order contribution of e M ( x ,x ^ — 1 is 

±ib&(x>,x",y) = I [(4 - 4')K - Wi) " Wi - - 2/2)] • (5.16) 

In first order in 6, (|5.14p gives: 

— ^2 s o(x,x';z)[i(x[~x / ()i(x 2 ~y 2 )-i(x' 2 -x 2 )i(x / (-y 1 )}h (x / ,x'')s (x / ',y;z), (5.17) 
x',x"eA 

where we have inserted some i's to bring it on the form (|5.12[) . Using (|5.17l) in expression (|5.14[) . 
we have that the linear contribution of 



Re ^T / dz /fd(^) V [S-hbQhb] (x,x) 

2\il\wuJo ±^ 



xen 

is: 



Re 4|^V' <1:/l " L ' ,: ' ' 



^ \r (z-uo)d 1 horo(z)[dih d 2 ro(z) - d 2 h diro(z)]d 2 ho {x,x), (5.18) 



xen 

1-00 4- Vi £1 i-» r\ Anon /-l fln ri*} r\\ +- It d An£ivo + to on/--l moo ^ It n otinn si- . 

dk. 



where we have suppressed the fc-dependence of the operators, and use the notation d\ = j^- and 

dl = m;- 
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5.1.1 All terms in fc-space 

Using the method above, we can calculate the traces of (I5.10[) , for a given z,z±lu in the resolvent 
set of Ho, explicitly, by simply inverting and differentiating known |Q| x |51|-matrices. 
Written in fc-space, these terms are: 

^2 [S-jibShb] (x,x) : 
xen 

( ^ k ^2 [( a i r o( z _ uo)d 2 dih - d 2 r (z - uj Q )d 2 h ^ r (z)d 2 h 
Jn ' xen 



2\n* 

+ [dir Q (z - uj )dih d 2 r (z) - d 2 r (z - Ld )dih dir a (z)) d 2 h 
+ r (z - ljq) (d 2 h d 2 r (z) - d 2 <9ih,o<9ir (.z)) d 2 h (x,x). 

l S -hbQhb] (x,x) : 



xen 
ib 



2\n* 



— d 2 fc \r {z - u )dih r (z) (dih d 2 r (z) - d 2 h Q dir {z)^ d 2 h (x,x) 
' xen 



53 [Q-5ibSj 2b } (x,x) 
xen 
ib 



ojT^T / ^ k ^2 [ r °( z _ (dih d 2 r (z - oj q ) - d 2 h Q dir (z - w )) dih r (z)d 2 h Q (x,x) 
' ' xen 

5.1.2 Collecting the terms 

Using that |f2||fT| = An 2 and inserting everyting into formula (|5.8j) we have: 



7 (i) = 

21 167T 3 W 



\ — / d 2 fcRe <p dz «/fd(-z) 
r J w J n . J v tf r 



xen 

2 



(5.19) 



^dir a (z - w )d 2 dih - d 2 r (z - wo^i^oj r (z)d 2 h 

+ (diro(z - uj )dih d 2 r (z) - d 2 r (z - Lo )dih dir (z)^ d 2 h 

+ r (z - uj ) (d 2 h d 2 r (z) - d 2 dih dir (z^j d 2 h 

+ r Q (z - u Q )dih r Q (z) (dihod 2 r (z) ~ d 2 h a dir a (z)j d 2 h 

+ r (z - loq) (dih d 2 ro(z - ui ) - d 2 h dir (z - w )J dih r Q (z)d 2 h 
+ {z ->• z + lu }] {x,x). 
This formula only contains known matrices and their derivatives. 

5.2 Consequences of the symmetry 

Now we want to prove (|2.4[) . The following lemmas will help us prove that all even Taylor 
coefficients of cr 21 (&) vanish. 

Lemma 5.2. For any n + 1-tuple of sites (x, x^\ . . . , x^) in A, it holds that 

n—l n—1 
tp(x,x<U) + V(x (m) .x^+V) + ip^^x) = J2 R ( x > x(m) > xim+1) )- ( 5 - 20 ) 
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Proof. By telescoping, the left hand side of the above equation can be rewritten as 

n-l 



[fix, x {m) ) + p(x ( - m \x^ m+1 '>) + ip{x^ m+1 \x) 



Identity (|1.10[) and the anti-symmetry of ip give (|5.20|) , 



□ 



Lemma 5.3. Given any (n+ 1) -tuple of sites (x, x^ , . . . , x^) in A, and an index 1 < r < n— 1, 
it holds that fi(x, x^ r \ x^ r+1 ^) is given by 



■r-l 



(x^l - X^ 



x (x - xW) + J2 ( x(r) ~ x [r+1) ) x (x (m) - x {m+1 ^) 



(5.21) 



Proof. A telescoping argument gives that 

fl(x, i< r ), x (r+1) ) = - [(a; (r) - x< r+1) ) x (x - x< r >)l 

(xW - x< r+1 )) X LxW + ^( 



which proves the lemma. 



□ 



Suppose that we want to detemine the n'th Taylor coefficient (n > 2) given like (|2.4j) . For 
z G p(Hb) put 5 = z — loq. The problem is basically to identify the coefficient to b n in matrix 
elements of the type 



' n 



n 



^ S b (z) (K b (z)) r j 1)b ^ S b (z) (K b (z)) r j 2)b 

\r=0 J \r=0 / 



(x,x), 



(5.22) 



as we did for n = 1 in expression (|5.8[) . Expression (|5.22[) can be expanded into a finite sum 
of terms all in the form [S b (2) (K b (z)) 91 ji, b S b (z) (Kb^)) 92 J2,b] (x,x) for some < <7i,<72 < 
Expanding one such term as a sum over products of integral kernels we obtain (to shorten notation, 
we write from now on x m instead of x' m )): 

[S b (z)(K h (z)) qi h,bS h {z)(K h {z)) q2 h,b](x,x) = ]T F 1 (x,x 1 )J' 2 (x 1 > x 2 )...i! , g (x9 > x) J 

X'eA,i=l...g 

where each •) is an operator kernel of either an S, a K or a j operator. The &-dependence of 
such an expression is always given in the form 



M(b) = e lb Mx,x 1 )+...+i P (x\x l+1 )+...+<p(x"^x)] M(b), 



(5.23) 



where M(b) is given by a convolution of kernels at b = 0, together with factors of the type 
^ e ?bti,(y,y ,y ) _ w here y, y" and y'" are consecutive convolution variables. Together with 

lemmas 15.21 and 15.31 it follows that the phases can be expanded in such a way that the only type 
of factors which can appear are of the form 

a(x,x',x",x"') = [(x - x') x (x" - x'")] z = (a* - x[)(4 - a£") - (x'{ - <)(a a - 4), 

where x and x' (respectively x" and x'") are consecutive variables in the convolutions. 

The coefficient of b n will consist of a finite number of convolutions, each of which having n such 
a factors. They will generate n factors of type j/2 — J/21 where y and y' are consecutive convolution 
variables. 
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We have to keep in mind that a 21 contains from the beginning an X2 — x' 2 coming from ] 2 . 
Thus the convolutions giving b n will all have exactly n + 1 factors like x 2 — x' 2 , where x and x' 
are consecutive convolution variables. Switching to the k space, these factors will be transformed 
into n + 1 partial derivatives with respect to k 2 . Remember that all matrix elements of /iq (k) and 
7q {z, k) are even functions of k 2 . Distributing 2p + 1 derivatives with respect to among these 
matrix elements will generate a global odd function in k 2 . When integrating with respect to ki 
over the symmetric first Brillouin zone, we get zero. Thus (|2.4[) is proved. 

6 Conclusions 

1. We constructed the conductivity tensor going through the Kubo-Greenwood formalism, pay- 
ing attention to the thermodynamic and adiabatic limits. Most physical papers completely 
ignore these issues, but this do not make them less important. Our proof of the thermody- 
namic limit is based on a simplified version of the geometric perturbation theory as developed 
in |13) for the Schrodinger case and then further developed in [14] . 

2. The magnetic perturbation theory cannot be avoided if one wants to control the growth at 
infinity induced by the constant magnetic field. Moreover, it provides us with a systematic 
method of computing derivatives of any order at 6 = 0. 

3. Remember that the Faraday rotation $ is proportional with <7 2 i(0) (denoted by a^i m 
formula (jl.2p . which gives the Verdet constant). In (15.191) we obtain a closed formula for a 2 ^ 
which will be the starting point of a further analysis of the dependency of the Verdet constant 
on temperature, chemical potential, density and spectral structure of a given material. For 
instance, the graphene is very interesting because at zero temperature the chemical potential 
lies exactly where the valence and conduction energy bands touch each other (see points K 
and K' in figure 2]) . The eventual lack of regularity of the Fermi surface can make the zero 
temperature limit nontrivial. 

4. The expression giving has an analytic extension in ujq to the whole complex plane except 
in zero and those real values for which the sets er(Ho) ± have common points with <t(Hq). 
But for certain particular models one can further extend the permitted regions of luq. In 
fact, it would be very interesting to study how a 2 ^ behaves when ujq comes close to those 
resonant values which induce transitions between different Bloch-bands. 

5. In the case when the magnetic field b generates a rational flux through the unit cell of Ho, 
then the spectrum of is absolutely continuous and consists of bands, but the elementary 
cells of H& become larger and larger when b becomes smaller and smaller. Nevertheless, one 
can compute 021 (&) in terms of the ^-dependent Bloch structure, see e.g. [15] . It would be 
interesting to compare this approach with our results. 

6. An open problem: take our graphene Hamiltonian H[f whose kernel is given in (jl.7p . and 
put a weak magnetic field on it through a Peierls phase. What happens with the spectrum 
of around the crossing of the valence and conduction bands, represented in figure 2]? 
Physicists claim that in that energy region the dynamics is close to the one generated by 
some zero-mass Dirac operator, and when we add a magnetic field it should create gaps. 

7. We note that the first term of (|2.ip dissapears only after the adiabatic limit (57 = 0). For 
graphene, unlike the usual Schrodinger operators, the commutator [j2,b,N, X^n] is not zero. 
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